function [result16]=wavProcess(filename)
f=wavread(filename,'native');
f=f(1:20480,1);
len=size(f);
len=len(1);
%inputData[i] - Alpha * (inputData[i - 1]);
pre(1)=f(1);
pre(2:len) =f(2:len)-0.95*f(1:len-1);
hm=hamming(1024);
hamwin=zeros(20,1024);
fftResult = zeros(20,1024);
magResult = zeros(20,1024);
logResult = zeros(20,1024);
IfftResult = zeros(20,1024);
for(i=0:19)
    hamwin(i+1,:)=hm'.*double(pre(i*1024+1:i*1024+1024));
end;

for(i=1:20)
    fftResult(i,:)=fft( hamwin(i,:));
end;
for(i=1:20)
    magResult(i,:)=( ( abs( fftResult(i,:) ) ) );
end;
for(i=1:20)
    logResult(i,:) =log10( magResult(i,:));
end;
for(i=1:20)
    IfftResult(i,:) =ifft( logResult(i,:));
end;
result=IfftResult';
result=result';
result=result(:);

result16 =zeros(20,16);
for(i=1:20)
    result16(i,:) =IfftResult( i,2:17);
end;
result16=result16';
result16=result16(:);



